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Thomas-Fermi density distribution in phase space representation, the total energy 
is described by a function of deformation parameters in both real and momentum 
space. Various thermodynamic quantities of a uniform dipolar Fermi gas are derived, 
and then instability of this system is discussed. For a trapped dipolar Fermi gas, the 
collective oscillation frequencies are derived with the energy- weighted sum rule method. 
The frequencies for the monopole and quadrupole modes are calculated, and softening 
against collapse is shown as the dipolar strength approaches the critical value. Finally, 
we investigate the effects of the dipolar interaction on the expansion dynamics of the 
Fermi gas and show how the dipolar effects manifest in an expanded cloud. 
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1. Introduction 

In recent years, atomic quantum dipolar gases have received much interest, for the 
simple reason that the anisotropic and long-range nature of the dipole-dipole interaction 
gives rise to a rich spectrum of novel properties to such systems. The theoretical 
study of dipolar Bose-Einstein condensates started in 2000. Properties of ground 
state [TJ |2], collective oscillations [3] 0], topological defects such as spin textures and 
vortex states [SI E] are studied. Moreover, when confined in optical lattice potentials, 
various quantum phases, such as ferromagnetism [7], and supersolid state [SI [9], etc. are 
predicted. Theoretical studies of dipolar Fermi gas have been carried out for ground 
state [10], excitatations [UJ, BCS superfluidity [12] and rotating properties [13J. A 
recent review of dipolar quantum gases can be found in Ref. [T3] , 

In experiments, Bose-Einstein condensation of chromium atoms, which possess 
a magnetic dipole moment six times larger than that of alkali atoms, have been 
realized [IS], [16]. The effect of dipole-dipole interaction in 52 Cr condensate is observed in 
its expansion dynamics [19] . Besides chromium, heteronuclear molecules [20] [2H [221 1231 
l24l [25] and Rydberg atoms [26], [27] [28] are also expected to interact via strong dipole- 
dipole force due to their large electric dipole moment, and their experimental realization 
is under way in a number of groups. 

In Ref. [2H], three of us studied the ground state properties of a dipolar Fermi gas by 
employing a variational Wigner function based on the Thomas- Fermi density of identical 
fermions. We showed that the dipole-dipole interaction induces a deformation of the 
momentum space distribution, and identified that such deformation arises from the Fock 
exchange term, which had not been paid particular attention in previous studies. The 
purpose of this paper is to extend the work of ref. [29] , and investigate the collective 
excitations and expansion dynamics of the dipolar Fermi gas. We want to emphasize 
that, due to the Pauli exclusion principle, the energy scales of a fermionic system is 
much larger than those of a Bose condensate. Consequently, the dipolar effects in Fermi 
gas only becomes significant when the dipole moment is very large. Our calculations 
show that for heteronuclear molecules with typical electric dipole moment on the order 
of one Debye, dipolar effects can be easily detected. While dipolar effects are usualy 
negligible in atomic Fermi gases 0. 

The content of the paper is organized as follows. In the next section, we present 
the model Hamiltonian and the total energy of the one-component dipolar Fermi gas 
under Hartree-Fock approximation. In section [3] we derive the total energy function 
in a uniform system with a variational ansatz of Fermi surface and compute various 
thermodynamic quantities of the system. Here we show how the Fock exchange 
interaction leads to Fermi surface deformation as well as the instability of the system. 
In section H] we turn our attention to a trapped system and investigate various modes 
of collective excitations using the sum-rule method, and show the softening of the 

X As pointed out in Ref. [14] . the magnetic dipole moment of chromium is equivalent to an electric 
dipole moment of 0.056 Debye. 
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excitation frequency as the interaction strength is increases towards a critical value. 
In section [5j we study the expansion dynamics of an initially trapped Fermi gas and 
show how the expanded cloud bears the signature of the underlying dipolar interaction. 
Finally, a summary is presented in section [6j 



2. Total energy functional in phase space representation 

We consider a single component Fermi gas of atoms or molecules with dipole moment 
aligned along the axial axis of a cylindrical harmonic trap. The Hamiltonian of this 
system is described by 

& = E + \m {co*(x> + y}) + + £ V M (n - (1) 

i L J i>j 

where m is the mass of fermions, and u p and uj z are the oscillation frequencies along 
the radial and axial axes, respectively. The dipole-dipole interaction of the last term in 
Eq. (pQ) is described by Vdd(?) = d 2 (l — 3 cos 2 9)/r 3 , where 6 is the angle between r and 
the dipole moment d. 

In the Hartree-Fock approximation, the total energy derived from Hamiltonian (pQ) 
can be written as the sum of the kinetic, trapping potential, Hartree direct and Fock 
exchange energies 

E = Ekin + E ho + E d + E ex (2) 



E d = \ f*Tj&f. U V dd{r - r')/(r, k)/(r', k') (5) 

E„. = - \ Jd 3 Rjd^sJ^J0^V M (s)e'*-«>-f(R,kif(RXt , (6) 

where we have introduced the Wigner function /(r, k) defined by the following 
transformation: 



' V ' r,= (2^/ <i3<;eik ' ( "' >/ ( £ ^ : ' k )' (?) 

where the one-body density matrix n(r, r') = ^ a ^(r)^! 1 "') is defined in terms of a 
complete set of single-particle wave function {^(r)}. In Eq. ([6]), we have introduced 
the center of mass coordinate R = (r + r')/2 and relative coordinate s = r — r'. For the 
ground state, the summation over single-particle states a goes from the lowest one up 
to the Fermi energy. 

In our work, we do not calculate the Hartree-Fock energy represented by Eq. ([2]) 
in a fully self-consistent manner, which will be a quite complicated task. Instead, we 
adopt a much simpler semiclassical approach and calculate the total energy by employing 
a variational ansatz for the Wigner distribution function based on the Thomas-Fermi 
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approximation, which assumes that the local Fermi surface has the same form as in 
homogeneous case at each spatial point. The ground state is then obtained by optimizing 
the Wigner function that minimizes the total energy. The details of this calculation can 
be found in Ref. [29]. In the present paper, we will focus on the dynamical properties 
such as the low-lying collective excitations and the expansion dynamics of the ground 
state. 

3. Equilibrium properties of a homogeneous dipolar Fermi gas 

It is instructive to first consider a homogeneous system (u p — lo z — 0) in a large box of 
volume V(= J d 3 r) with number density rif, as this will provide important insights into 
the trapped system to be studied later. 

We introduce the following number-conserving variational ansatz for the Wigner 
function 



to the Fermi momentum. The parameter a characterizes the deformation of the Fermi 
surface: a > 1 (< 1) corresponds to an oblate (prolate) Fermi surface. The physical 
origin of the Fermi surface deformation can be attributed to the anisotropic nature of 
the dipolar interaction. 

Given the ansatz Eq. ([8]), the total energy of the homogeneous system can be derived 

as 



is the "deformation function" [29] and is illustrated in Figure [TJ I(x) decreases 
monotonically from 4 to —2 as x increases from to oo, and passes through zero at 
x = 1. The first and the second term in the square bracket of Eq. ([9]) represent the 
kinetic and the Fock exchange energy, respectively. For the homogeneous system, the 
Hartree direct energy vanishes. 






where C\ = 3(67r 2 ) 2 / 3 /10, C^d — md 2 n x j z /h 2 is the dimensionless dipolar interaction 
strength, and 
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Figure 1. Deformation function I(x) as a function of x. 
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Figure 2. Optimal deformation parameter ao as a function of Cdd- For a molecular 
Fermi gas with electrical dipole moment d — 1 Debye, molecular mass m = 100 a.m.u. 
and density rif = 10 13 cm~ 3 , we have Cdd ~ 3.2. 

Under this variational approach, the ground state is determined by the stationary- 
condition for the total energy of Eq. ([9]) with respect to parameter a: [de /da] a= = 0. 
The optimal value ao is shown in Figure Has a function of the dipolar strength Cdd- 
For free fermion systems, momentum density distribution is spherical, i.e., «o = 1 at 
Cdd = 0. As the interaction strength increases, «o decreases, which means that the 
momentum density distribution becomes more prolate in shape. In other words, the 
Fermi surface is stretched along the direction of the dipoles. 

Once we have the energy of the system as represented by Eq. (Q, we can easily 
obtain other important thermodynamic quantities. Here we provide our calculation for 
the pressure P, compressibility K and chemical potential /i: 
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These quantities are illustrated in Figure [31 One can see that P, 1/K and \x all 
monotonically decrease as the dipolar interaction strength increases. In particular, when 
the inverse compressibility (i.e., the bulk modulus) becomes negative, the system is no 
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Figure 3. Chemical potential y, pressure P, and inverse compressibility or 
bulk modulus l/K as functions of Cdd- All quantities are normalized to their 
corresponding values in the non-interacting limit: yo — (5Ci/3)h 2 n 2 ^ 3 jra — Ef, 

P q = (2C 1 /3)h 2 n 5 f /3 /m, and 1/K = (lOd/tytfn 5 / 3 /m. The vertical line indicates 
the critical dipolar strength beyond which the system becomes unstable against 
collapse. 

longer stable against collapse. Our calculation indicates that the critical dipolar strength 
is about Cdd = 3.23. 

4. Collective oscillations of trapped dipolar Fermi gas 

Let us now turn our attention to the trapped dipolar Fermi gas. 

First, to obtain the the total energy of Eq. ([2]), we introduce the following ansatz 
for the Wigner function: 

/(r , k ) = e(4-f-^-|(^ + 3), do) 

where p 2 = x 2 + y 2 , and a^ = y/h/ (mu) with u = (cu 2 ^) 1 / 3 . The variables j3 and 
A represent the deformation and compression of the spatial density distribution of the 
system, respectively. When we take a — 1, /3 — (uj p /uj z ) 2 ^ 3 , and A = 1, this trial function 
is consistent with the Thomas-Fermi density of a free Fermi gas in the harmonic trap. 
Fermi wave number kp is related to the number of fermions as 

N = j d 3 rn(r) = J £r J H^r.k) = !« (11) 

Substituting Eq. ffTUj) into Eqs. (j3J), (jl]), (j3J), and (jSJ), we obtain the total energy in 
units of N^hw as [29] 

e(a, (3, 7) = Ni/3huj = £kin(a, A) + e ho {(3, A) + e d (f3, A) + e ex (a, A) (12) 

e ho = — ( — + ^2 ( 13 ) 



efein = CiA (^2a + -jjj 



a V P Ph 

e d = N 1 ^c d dC2l(P)X 3/2 , e ex = -iV 1 / 6 c dd c 2 J(a)A 3 / 2 (14) 

where (3 = (cu p /uj z ) 2 ^ 3 represents the trap aspect ratio, c\ = 3 1 / 3 /2 8 / 3 ~ 0.2271, c 2 = 
2 10 /(3 7//2 -5-77r 2 ) ~ 0.0634, and c dd = d 2 /(huia^) is the dimensionless dipolar interaction 
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Figure 4. (a) Stability phase diagram in the space of the trap aspect ratio and 
the dipolar interaction strength, (b) Different energy terms, in units of N 4 / 3 huj as 
functions of the trap aspect ratio (3q at N 1 l 6 Cdd = 1.5. 



strength for the trapped system. The momentum space deformation parameter a, as in 
the homogeneous case, appears only in the kinetic and the exchange energy terms, both 
of which are independent of the spatial deformation parameter (5. This indicates that 
the momentum space distribution of the trapped system will also be elongated along 
the direction of the dipoles, regardless the geometry of the trapping potential. On the 
other hand, (3 appears only in the potential energy and the Hartree direct energy terms. 

The ground state is determined by the stationary condition for Eq. (|T2l) with respect 
to the three variables a, (3 and A: de/da = de/df3 = de/dX = 0. From the last condition, 
we can see that the energies of the dipolar Fermi gas satisfy the Virial theorem: 

2efcj n — 2eft + 3(e<2 + e ex ) = 

In addition, the ground state has to satisfy the stability condition: The energy surface 
in the coordinates (ot,/3, A) has to be a convex downward at the stationary point. If no 
values of (a, /3, A) can be found to satisfy both the stationary and the stability conditions, 
the systems is considered to be unstable against collapse [29]. This procedure leads to 
the stability phase diagram as shown in Figure H](a). Just as in the homogeneous 
case, the trapped dipolar gas is only stable for dipolar interaction strength below a 
critical value. In Figure IU(b), we show the different energy terms [Eqs. (1131) and (TT4j) ] 
as functions of (3q at N l l % Cdd = 1.5. Several features are worth pointing out: (1) The 
exchange energy is always negative, as in the homogeneous case, regardless of the trap 
geometry; whereas the sign of the direct energy depends on trap geometry: > for 
(3q < 1 (oblate trap) and < for /3 > 1 (prolate trap). (2) Both the kinetic and the 
trapping energies depend on trap aspect ratio. By contrast, for non-interacting system, 
when expressed in the same units, we have = €ho = 3ci ~ 0.68 independent of /So- 
Next, we derive the collective oscillation frequency for several low- lying excitation 
modes of the system using the sum rules in the present formulation [301 El]- In this 
approach, we represent the excitation frequency Q using the first and third energy- 
weighted moments of the strength function for a given transition operator F: 
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Si = - £o)|(^|0>| 2 = i<0|[F, [H, F]]\0) (16) 

S 3 = Y,t E » - ^o) 3 |(H^|0)| 2 = -\(0\[[H,F}, [H, [H,F]]]\0), (17) 

where \v) denotes the z/-th eigenstate of the Hamiltonian with eigenenergy E v . 
For our purpose, we choose the one-body operator as 

N N 

where £ and £ are certain parameters. A collective oscillation is compressive when 
£ and ( have the same sign, and quadrupolar when they have opposite signs. The 
natural monopole and quadrupole operators correspond to £/£ = 1 and £/£ = —1/2, 
respectively. 

Using Eq. (ITS"]) , the collective excitation frequency Q in Eq. ([15]) in the present 
formulation can be shown to be 



n _ 4(6, op e + e^c 2 ) + ^ - oc + ^(4e - cxe - o + c{i - o s 



(19) 



1/ \ 2 A,<9ed <9e ex \ „ 2/ 9 <9 2 e d od 2 e ex \ , . 

.4=^+^), e =i\fw +c? -m (20) 

where = 2c\(3q/(\(3) and = CiP 2 /(\(3q) are the radial and axial components of 
the trapping energy, respectively [see Eq. ffl~3]) ]. The excitation frequency f2 is calculated 
by substituting the variational parameters (a, (3, A) at the stationary point of the total 
energy (fl2|) . 

From Eq. ( 1191) we can easily find the excitation frequencies of the monopole and 
quadrupole modes, which have the following expressions: 



n. tf _ u I Ae ho + 3A pi) 
(-hop I Po + (3q 



n / 4e hop + 16e fto2 - 12.4 + 623 + 9C 

"° " " V + • (22) 

The corresponding frequencies for non-interacting system can be recovered from 



Eqs. (EU) and ([22]) as Q M0 = ^ 12u 2 u 2 / {u 2 + 2u 2 ) and Q Q0 = ^12^/(2^ + w 2). 

Figure [5] shows the excitation frequency of the monopole mode Qm and the 
quadrupole mode VLq. As can be seen in Figure IU(b), the total interaction energy 
(e ex + €d) is positive (in other words, the overall dipolar interaction is repulsive) for 
oblate traps (/3 < 1) which makes the atomic cloud more compressible, hence Vt M 
is increased compared to its non-interacting values. For prolate traps (/3q > 1), the 
opposite will be true. This is consistent with the result shown in Figure Ufa). The 
quadrupole mode frequency Qq, on the other hand, exhibits a roughly opposite trend. 

To account for the hybridization of different modes, we parameterize £ and ( in 
Eq. (TT81) as £ = sin# and ( = cosd, with < 6 < it. We then investigate the minimum 
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Figure 5. Excitation frequencies of the monopole mode Q,m (£ = C = 1) an( i the 
quadrupole mode Qq (£ = 1, £ — — 2) as a function of /3q. The frequencies are 
normalized to the corresponding values of the non-interacting system. 
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Figure 6. The minimum excitation frequency f2 m i n (a) and the angle that minimizes 
£1(6) (b) as functions of the interaction strength N 1 / & Cdd up to the critical values 
against instability. f2 m in is normalized to f2o, the corresponding frequency for the 
non-interacting system at each f3 : fl /uj = 1.7, 2.0 and 0.2 for /3 = 0.8, 1.0 and 10, 
respectively. The critical values are A^ 1/6 c drf = 2.433, 2.166 and 1.603 for f3 = 0.8, 1.0 
and 10, respectively, see Figure |U 



value of the excitation frequency Q(6) given by Eq. ffl9|) . The collective oscillation will 
be dominated by the compression mode for < 9 < tv/2 and by the quadrupolar mode 
for 7r/2 < 9 < 7i. Moreover, 9 = n/2 represents a radial mode, and 9 = an axial 
mode. The natural monopole and quadrupole operators correspond to 9 = 7r/4 and 
9 = 7r — arctan(l/2) w 0.857T, respectively. 

Figure HD^a) shows the minimum excitation frequency fi m m as a function of the 
interaction strength N 1 / & Cdd up to the critical value, while Figure E](b) shows the angle 
9 that minimizes Q(9). For the spherical trap with [3q = 1.0, the excitation frequency 
decreases monotonically as the interaction strength increases, and the minimum-energy 
mode is the monopole mode. For the prolate trap with [3q = 10.0, the minimum-energy 
mode is dominated by the axial mode with 9 m as the axial axis represents the direction 
of the soft confinement. Similarly, for the oblate trap with (3q = 0.8, the minimum- 
energy mode is dominated by the radial mode with 9 ~ tt/2 as the radial direction now 
becomes the soft axis. However, as the interaction strength increases towards the critical 
value, in both of these cases, the minimum-energy mode shifts towards the monopole 
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mode, and we clearly see the tendency of the softening of the collective mode, indicating 
the approaching of the collapse instability. We note that, in particular for the case of 
flo = 0.8, fimin does not completely decrease to zero at the critical value. This could be 
due to the calculation of the average frequency of the collective oscillation by sum-rule 
method. Deeper insights into collective excitations may be obtained from microscopic 
approaches such as the random-phase approximation [32j [33] . 



5. Expansion dynamics 



We now turn to the expansion dynamics of an initially trapped dipolar Fermi gas. This 
study is important as in most cold atom experiments, the atomic cloud is imaged after a 
period of free expansion. Furthermore, the expansion dynamics may bear the signature 
of the underlying interaction. The dipolar effects in chromium condensate are first 
observed in the expansion dynamics p2J [TFl [18] . 

Our starting point is the Boltzman-Vlasov equation: 

+ — + T^rU(r,k,t)) ■ T-f (r,k,t) 



dt \m hdk ) dr 

-~tf(r,M).;|-/(r,M) = J (23) 
n or ok 

where the effective potential U includes both the external harmonic trap potential Uho 
and the mean-field potential due to the dipole-dipole interaction: 

dV n(r', t) V dd (r - r') - J j^y 3 V dd (k - k') / (r, k', t) , (24) 

where V dd (k) = (And 2 /3){?>k 2 z /k 2 — 1) is the Fourier transform of V dd {v). Note that the 
k-dependence of the effective potential U originates exclusively from the contribution 
of the exchange interaction, i.e., the last term at the r.h.s. of Eq. ( |24l . 

To study the dynamics, we shall make use of the scaling transformation 

/(r,M) = /o(R(t),K(t)), 

Ri{t) = n/bi(t) , Ki(t) = b i (t)k i -mb i (t)r i /h, (i = x,y,z) 

where /o represents the equilibrium Wigner distribution function obtained in previous 
section, whose form is given by Eq. ffTUj) . and 6« are the dimensionless scaling parameters. 
This scaling approach has been used previously to study the expansion of Fermi gases 
[2U [35] and Bose- Fermi mixtures [36J. 

From the Boltzman-Vlasov equation, we can derive the equations governing the 
scaling parameters [SI [36] 

, 2 



7}(b,b) 7Kl,0) 



0, (25) 



with 7j = ujj/uo, e dd = N l l & c dd and (R 2 ) = J d 3 R R 2 no(H) with no being the equilibrium 
density. The second and third terms in Eq. (1231) represent, respectively, the restoring 
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force and the kinetic energy. Collecting all contributions from interaction we have 

dn (R') 



7~(b, b) = / d 3 Rd 3 R' RjW{b; R - R')n (R) 



dR> 



+ / d 3 Rd 3 Kd 3 K' RjKj W (b; K — K') 



E 



g/o(R,K) 



9/o(R,K0 bAdf (R,K') 



dRi 



2n 



dK' 



+ J d 3 Rd 3 Kd 3 K'R j K j f (R, K') 



E 



(b; K - K') 



0/o (R,K) , 6i6 f a/ (R,K) 



2tt 



dKi 



(26) 



where W(b; R) 



(tfx 2 +&Y 2 +tfz*) 5 / 2 is the dipole-dipole interaction potential under the 
scaling transformation and W(b; K) represents its Fourier transform. Given the Wigner 
function in Eq. (ITU]) , we obviously have b x = b y = b p as the free expansion will not change 
the cylindrical symmetry. Moreover, the integrations for terms involving b{ in Eq. (1261) 
vanish, so that Tj reduces to a function of b only. The analytical expressions for Tj (b) 
can be found as 



T p (b) = qb-% 1 [d p (p- 3 %/b z ) - d p (a- 3 %/b z )} , 
%{b) = qbfb- 1 [d z (p- 3 %/b z ) - d z {a- 3 ' 2 b p /b z )] , 

where q = 1024(3A 3 ) 1 / 2 /(2835vr 2 ) and the functions dj are defined as 

d p (x) = (1 - x 2 Y 2 [2 - 7x 2 - 4x 4 + 9x 4 g(x 



(27) 



d z (x) = 2(1 - x 2 Y 2 [1 + 10x 2 - 2x 4 - 9x 2 g(x) 

with g(x) = tanh -1 -^/! — x 2 /\/l — x 2 . We note that dj are all monotonically decreasing 
functions of x and bounded between 2 and —4 for x G [0, 00). 

Here we focus on the time evolution of the atomic cloud aspect ratios in real and 
momentum spaces which are defined, respectively, as 



K r (t) 



(r P (t) 2 ) 



and K p (t) 



(P P (t) 2 ) 
(P z (t) 2 ) 



(r z (t) 2 ) 

where initially the system is prepared in the ground state inside the external trap. 
Straightforward calculations yield that 

Kr{t)=/3- 3/ X(t)/b g {t) 



Kp(t) 



P\ 2 a% 2 + a 2 b 2 



1 1/2 



(3X 2 b z 2 + a 2 p 3 b 2 z 

The initial cloud aspect ratios are determined by the ground state Wigner function 
and can be easily shown to be n r (0) = (3~ 3 ^ 2 and k p (0) = a 3 / 2 . To study the expansion 
dynamics, we turn off the trapping potential at t — and the cloud starts to expand. We 
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Figure 7. Cloud aspect ratio during time of flight in both momentum space (dashed 
lines) and real space (solid lines) for 0o = 0.5 (a), 1 (b), and 5 (c). In each figure, the 
dipolar interaction strength are N^^Cdd = 0.5, 1.0, and 1.5, in descending order. 

then solve for b p (t) and b z (t) using Eq. (1251) with the restoring force term jjbj removed 
and with the initial conditions b p (0) = b z (0) = 1. Before presenting our results, we recall 
that when the exchange interaction is ignored, the direct dipolar interaction always tends 
to stretch the cloud along the direction of dipole moments in both real and momentum 
spaces [35] . 

Figure [7] displays several examples of the cloud aspect ratio during time of flight for 
different trap geometries. As expected, asymptotically, the aspect ratios in momentum 
and real spaces become equal to each other, i.e., K r (oo) = k p (oo) = k^. A notable 
feature is that, regardless of the initial trap geometry, the shape of the expanding cloud 
eventually becomes prolate as < 1. This result is in stark contrast to the expansion 
dynamics of a dipolar condensate whose asymptotic aspect ratio is sensitive to the initial 
trap geometry [T6l [TTl [18] . Furthermore, the interaction effects during the time of flight 
is also evident from Figure [7J Had interaction been ignored, the expansion would have 
become ballistic with k p a constant in time. Figure [7(b) indicates that the expansion 
is essentially ballistic for an initial spherical trapping potential, as for such traps, the 
interaction energy is rather weak as shown in Fig. H^b). 

That the expanded cloud eventually becomes prolate in shape is also obtained in 
Ref. [35] when the exchange dipolar interaction is ignored, which indicates that the 
effect of the exchange interaction during the expansion is not very important. This 
is consistent with Fig. H^b) which shows that, except for nearly spherical traps, the 
magnitude of the direct energy is in general much larger than that of the exchange 
energy. However, we want to emphasize that the exchange term is crucial for the 
equailibrium momentum distribution inside the trap: Without the exchange term, the 
momentum distribution would be isotropic for any trap geometry. To get a closer look, 
we compare in Fig. [8] with the initial momentum space aspect ratio k p (0) which 
characterizes the momentum distribution for the ground state in the trap. The initial 
momentum distribution is always prolate in shape as k p (0) < 1. In general, the effect 
of the interaction during the expansion, with the dominant contribution from the direct 
term, is to further enhance this anisotropy such that < k p (Q). Exceptions may occur 
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Figure 8. (a) The dipolar interaction strength dependences of asymptotic aspect 
ratios Kqo (solid lines) and the initial momentum space aspect ratio k p (0) (dashed 
lines) for various trap aspect ratio /3q's. (b) The difference between the asymptotic 
aspect ratio and the initial momentum space aspect ratio for /?o = 1. 

for nearly spherical traps, for which one may have > /%,(0) as shown in Fig. E^b). 
However, this effect is very small since, as we have already mentioned earlier, the total 
dipolar interaction is weak for such traps. 

6. Summary 

In summary, we have studied the properties of dipolar Fermi gases in both homogeneous 
system and in the cylindrical harmonic trap with the dipole moments oriented along 
the symmetry axis. The total energy functional of this system is derived under the 
Hartree-Fock approximation. The one-body density matrix in the energy functional is 
obtained from a variational ansatz based on the Thomas-Fermi density distribution in 
phase-space representation, which accounts for the interaction-induced deformation in 
both real and momentum space. Our calculations show that deformation of the spatial 
density distribution comes from the Hartree direct energy term, while deformation of 
the momentum density distribution arises from the Fock exchange energy term. Note 
that the exchange term, a consequence of the anti-symmetry of the many-body fermionic 
wave function, does not appear in Bose-Einstein condensate. 

We have calculated several thermodynamic quantities such as the pressure, the 
compressibility and the chemical potential of the homogeneous system and investigated 
the low-lying collective excitations of a trapped dipolar Fermi gas using the sum rule 
method for various trap geometry and interaction strengths. We observe the softening 
of the collective excitations as the interaction strength approaches the critical value for 
collapse. 

Finally, we have studied the expansion dynamics of the initially trapped system. 
We show that, in stark contrast to dipolar condensate [T6], [T71 ITB] . the atomic Fermi 
gas will eventually become elongated along the direction of the dipoles regardless of the 
initial trap geometry. This feature makes it convenient to detect the dipolar effects in 
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Fermi gases. 
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